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We address an old problem in lattice gauge theory - the computation of the spectrum and wave functions 
of excited states. Our method is based on the Hamiltonian formulation of lattice gauge theory. As strategy, 
we propose to construct a stochastic basis of Bargmann link states, drawn from a physical probability density 
distribution. Then we compute transition amplitudes between stochastic basis states. From a matrix of transition 
elements we extract energy spectra and wave functions. We apply this method to U( 1)2+1 lattice gauge theory. 
We test the method by computing the energy spectrum, wave functions and thermodynamic al functions of the 
electric Hamiltonian of this theory and compare them with analytical results. We observe a reasonable scaling 
of energies and wave functions in the variable of time. We also present first results on a small lattice for the full 
Hamiltonian including the magnetic term. 

PACS numbers: 



I. INTRODUCTION 

Since the invention of lattice gauge theory (LGT) by 
Wilson in 1974 UJ], it has been customary to compute 
Euclidean 2-point functions (and likewise n-point func- 
tions), which is the expectation value of two field oper- 
ators in the vacuum-to-vacuum transition amplitude, 

(n,f = +^\'¥{x)'¥{y)\Q.,t = -00) . (1) 

The term Euclidean means that one works in imaginary 
time t —i' —it. It is important to note that the ground 
state wave function Q. is actually not known in LGT. 
The ground state wave function emerges in the above 
amplitude via the Feynman-Kac limit through projec- 
tion from an arbitrary state |E) onto the vacuum state in 
the limit of large Euclidean time, 

exp[-:// t/h] |S) - lim exp[-E,,,.t/ri] \Q.) {Q.\Z) . (2) 

If one wants to construct the wave functions of the 
ground state or excited states, one needs information 
beyond the vacuum-to-vacuum amplitude. Thus let us 
consider transition amplitudes 

(Tv,f = r|T^,f = 0) = (Tv|exp[-iWr/ri]|T^) (3) 

between initial and final states taken from some set 
of states |Tv),V = 1,...,A^. Here |Tv) denotes time- 
independent Bargmann-link states, that is, a configu- 
ration of link variables U assigned to each link ij on 
the whole spatial lattice. As a reminder, in the case 
of a chain of coupled oscillators, a Bargmann state de- 
notes the ensemble of displacements of the particles 
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from their resting positions. Hence Bargmann states can 
be interpreted as the analogue of position states in quan- 
tum mechanics. 

It is crucial to choose the states |Tv) to be physi- 
cally relevant and important. In the history of quan- 
tum mechanics it has been a tradition to use basis states 
obeying mathematical properties, like being orthogonal 
functions and forming a complete basis. Examples are 
Fourier functions, Hermite functions, etc. These func- 
tions have no relation to the particular physical system, 
e.g. a hydrogen atom. This is a luxury which one can 
afford when solving a system with few degrees of free- 
dom. However, in many-body physics and quantum 
field theory, with a huge (infinite) number of degrees 
of freedom, one better chooses a basis suited to the par- 
ticular physical system. 

There is a long tradition in many-body nuclear 
physics, condensed matter physics and elementary par- 
ticle physics aiming to compute wave functions and the 
energy spectrum from matrix elements of a Hamilton 
operator in a suitable basis. This approach has been suc- 
cessful in special cases, e.g., when perturbation theory 
is applicable (nuclear shell model, Kondo effect |01), 
or in situations where few degrees of freedom describe 
physics (Schwinger model in 1h-1 dimensions, renor- 
malization group Hamiltonian for critical phenomena). 
However, in general, this approach resulted in mitigated 
success. The problem is due to the choice of basis 
states. This can be understood by the example of diago- 
nalization of a real symmetiic matrix of large but finite 
rank, which is ill-conditioned (large difference in order 
of magnitude between largest and smallest eigenvalue). 
Any set of orthogonal basis states will yield a few lead- 
ing order eigenvalues. However, only basis states close 
to the eigen basis will resolve eigenvalues beyond lead- 
ing order. 

Carrying over this lesson to many-body physics and 
field theory, it is not surprising that perturbatively con- 
structed basis states (e.g. Fock states) in general are not 
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adequate, i.e., those states do not reflect the physicafly 
important degrees of freedom. This goes hand in hand 
with another potentially serious problem: Is there a con- 
trollable cut-off? The Tamm-Dancoff cut-off gives no 
estimate of the remainder Also, it can happen (e.g., in 
the coupled cluster method) that one has a quite large 
set of basis states (in the order of 10^ to 10^) yielding 
many tiny matrix elements all of about the same order 
of magnitude. Where to draw the line? 

We suggest a solution to such problems lying in a 
combination of two strategies: (i) Use stochastic tech- 
niques to sample states from a huge variety of possibil- 
ities, (ii) Use information from the physical system to 
guide the sampling. The above principles (i,ii) mean to 
do Monte Carlo with importance sampling. This pro- 
cedure has proven most successful in LGT to compute 
path integrals via generation of equilibrium path con- 
figurations. The stochastic basis states will be closely 
related to those equilibrium path configurations. 

From computation of transition amplitudes between 
stochastic basis states, one obtains a spectrum and wave 
functions of an effective Hamiltonian - the so-called 
Monte Carlo Hamiltonian - being valid in a low en- 
ergy, respectively low-temperature window. The idea 
of the Monte Carlo Hamiltonian has been suggested in 
1999 in Ref. lUt]. Its working in quantum mechanics 
has been demonstrated in a number of cases. For the 1- 
D harmonic oscillator, energy spectrum, wave functions 
and thermodynamical functions have been found to be 
in good agreement with the exact results IH-jll. Similar 
results have been obtained for uncoupled as well as cou- 
pled harmonic oscillators in 2-D 10131 and 3-D Jst. This 
has been extended to a variety of other potentials in 1-D 
like y oc ^,Voc \x\/2, and V Q{x)x HI, as 

well as the l/r Coulomb potential with a singularity at 
the origin yjj. The Monte Carlo Hamiltonian has been 
applied to the Yukawa potential V = —Vq exp{—ar)/r 
in the search for a critical value of above which 
no bound states exist lfl2ll . In field theory, the Monte 
Carlo Hamiltonian has been applied to the ( 1 + 1 ) Klein- 
Gordon model for the computation of the spectrum and 
thermodynamical functions il3Hl6ll , and likewise to the 
( 1 + 1 ) scalar model for the computation of the spectrum 
and thermodynamical functions fl7l - [l9ll . A first step 
towards the Monte Carlo Hamiltonian in lattice gauge 
theory has been made in ll20ll by computing transition 
amplitudes of U(l) gauge theory. 

In lattice gauge theory the U(l) model in (2 + 1) 
dimensions has been investigated since the early days 
using Euclidean Monte Carlo methods. Bhanot and 
Creutz D'Hoker IH, Ambjorn et al. Q com- 
puted the Wilson loop. Sterling and Greensite ll24ll 
computed the string tension from energy differences 
using external sources. The Hamiltonian formulation 
of the U( 1)2+1 model has been investigated using a 
variety of different methods. We can refer to the 
Quantum Monte Carlo methods lUsj-HH, the projector 



Monte Carlo method ll29l [30I1 . the ensemble projector 
Monte Carlo method Il3lll . the Green's function Monte 
Carlo method lf32ll . the Langevin technique ifssll and 
the guided random walk method ||34!l . Furthermore, the 
model has been studied via the t-expansion method Issl 
[36I1 . the block renormalization group rnethod iIstIi . the 
correlated basis function method |38||, the coupled- 
cluster expansion ll39l l40ll . and series expansions llill 
I42I1 . The U( 1)2+1 model has been also investigated in 
the finite lattice Hamiltonian approach by Allesandrini 
et al. ll43ll . and Irving et al. [44]. Recent results from the 
path integral Monte Carlo approach have been reported 
iniH. 

In this work we construct the Monte Carlo Hamil- 
tonian in U( 1)2+1 lattice gauge theory and apply it to 
compute the energy spectrum of excited states, the cor- 
responding wave functions and thermodynamical func- 
tions. The goal of this article is to show the working of 
the method applied to an Abelian, but non-trivial model 
in lattice gauge theory. We consider a few small spatial 
lattice volumes up to 8^. Because of the smallness of the 
lattices we do not consider here the quantum continuum 
limit (a 0). Rather, we try to give a careful analy- 
sis of the origin and size of errors, which determine the 
limitations of the method. We believe the potential of 
the method lies in its capacity to determine energies and 
wave functions for a number of excited states. 



II. MONTE CARLO HAMILTONIAN IN QUANTUM 
MECHANICS 

Let us consider in 1-D a massive particle moving in a 
potential V{x). Let S denote its classical action and let 
!}{ denote its quantum mechanical Hamilton operator. 
First we construct a stochastic basis, adapted to the par- 
ticular physical system. We do this by a random draw of 
nodes x, guided by the action S, respectively, Hamilto- 
nian operator !}{ of the system. We draw from the prob- 
ability distribution, given by the transition amplitude in 
imaginary time 

= Z {x\exp[-^H:T/n]\x = XMt) , (4) 

where Zl is a normalization factor such that / dx P{x) = 
1. The function P{x) is positive, because the tran- 
sition amplitude, expressed as path integral is given 
by a positive measure [dx] and a positive weight fac- 
tor exp[— S'euci/^] ('^Euci denotes the Euclidean action). 
For example, in case of the free Hamiltonian, choosing 
JCinit = 0, the distribution is a Gaussian, 

P{x) = (Aiexp[-:?4mrA]|xi„i, = 0) 

Here, i^in commutes with generators of translation and 
the measure dx is translationally invariant. In case of 



tion, / dxP{x) = 1, implies 
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FIG. 1: Quantum mechanics. Top panel: typical paths from 
path integral. Bottom panel: Scheme of box functions. The 
position of the center of any box function is drawn from a 
Gaussian distribution. They form a stochastic basis. 



the harmonic oscillator, P{x) is also given by a Gaus- 
sian. In order to draw samples from the distribution P{x) 
given by Eq. (|4|, one can use the path integral in imag- 
inary time. For example. Fig. [I](top panel) schemati- 
cally shows the sampling of points ,y, from the distribu- 
tion P{x) with the time parameter t ~ T, via the path 
integral going from f = to f = T, by intersecting paths 
at f = r/2. In such a way one obtains a set of nodes 
{xi,. . . ,xn}- This yields a stochastic basis of (non- 
normalizable) position states . . . , \xn). In order to 
eventually compute normalized wave functions of the 
Hamiltonian, one needs a basis of normalized Hilbert 
states. We introduce a stochastic basis of square inte- 
grable Hilbert states ("box states"), defined by 



bi{x) = 



hi if Xi 
else 



- Axi/2 <x < Xi + Axi/2 



(6) 



Those states are "quasi-localized", i.e., the function 
bi{x) is located around the node x; and has a "small" 
width Axi. The normalization condition of the distribu- 
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Thus we define 



Ax,- = 



NPixi) 



(7) 



(8) 



The normalization condition of the box functions, 
/ dx bj{x) ~ 1, yields 



hfAxi = 1 . 



Thus one obtains 
hi 



= VNP{xi) 

/Axi 



(9) 



(10) 



Fig. [T] schematically shows the stochastic basis of box 
states. 

Next, let us consider the computation of transition 
amplitudes in imaginary time between a pair of initial 
and final position states \y) and |z), respectively. 



M,,,^{z\exp[-^T/h]\y) 
[dx]&x-p[-~S[x]/h] 



(11) 



v,0 



Here S[x] denotes the Euclidean action of a path x{t) 
going from y to z- 



S\x] 



T I 

dt —mx^ -\-V{x) 
2 



(12) 



The Monte Carlo method with importance sampling is 
suited and conventionally applied to estimate a ratio of 
path integrals. Thus we suggest to calculate the matrix 
elements M-,, by splitting the action 



S = Sq + Sv 



JO 2 
and to express M~j as 



dt V{x), (13) 



Mz..y{T) = f [dx] exp[-5oW/ri] 



>:0 



z.T 



J[dx] exp[-5v/[x]/ri] exp[-5o[x]/;i]|^,;o 



z.T 



(14) 



Here MiJ denotes the matrix elements corresponding to 
the free action ^o, which is known analytically, that is, 

exp\-^{z~yf]. (15) 
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The term R-,y{T) denotes the ratio of path integrals. 
Such ratio can be computed by standard Monte Carlo 
methods with importance sampling, by treating O = 
exp[—Sv/Ti] as observable. 

Next, we want to compute transition amplitudes in 
imaginary time between a pair of initial and final states 
taken from the stochastic basis of box states. Thus we 
build the finite, real, symmetric matrix 

Mji{T) - {bj\exp[-l^{T/rl]\bi) 

rXj+Axj/2 rXj+Axj/2 

= hjhi / dz / dyM,,,{T) , 

Jxj—Axj/2 Jxj—Axj/2 

(16) 

where j run over 1 , . . . , A^. Using the property that box 
states are quasi-local, and assuming the transition am- 
plitudes vary little in the domain covered by the box 
functions, we may approximate Mjj in the following 
way 

Mj,{T) « lijh,AxjAx,Mi%,iT)xR,^^,iT) 
^ ^/MjMi Mi%- {T) X R^j^- {T) . 

(17) 

After having obtained the matrix M(T), given by 
Eq. ( fTTl ). one can extract eigenvalues and wave func- 
tions, which define an effective Hamiltonian. M{T) is 
a positive and Hermitian matrix, under the assumption 
that !}{ is a Hermitian operator Elementary linear alge- 
bra implies that there is a unitary matrix U and a real, 
diagonal matrix CD such that 

MiT) =U^ 'D{T)U. (18) 

On the other hand, projecting 9{ onto the subspace 
generated by the first states of the basis \bi), and using 
the eigen representation of such Hamiltonian, one has 



Mji{T) = Y^{bj\Ef)e 

k=\ 



'Ef\bd, 



and we can identify 



eff\ 



(19) 



(20) 



Thus algebraic diagonalization of the matrix M{T) 
yields eigenvalues 'Dk{T), k = 1,---,N, which by 
Eq. ( |20l ) determines the spectrum of energies. 



(21) 



The corresponding k-th eigenvector can be identified 
with the k-th column of the matrix l/Jj^. Hence from 
Eq. ( l20b we also obtain the wave function of the k-th 
eigenstate expressed in terms of the basis \bi). Thus 
starting from the matrix elements Mij{T) we have ex- 
plicitly constructed an effective Hamiltonian in a diago- 
nal form 



effi 



(22) 



k=l 



III. U(l)2+i LATTICE GAUGE THEORY 

The strategy to construct an effective Hamiltonian, 
outlined above for 1-D quantum mechanics, can be 
translated to lattice gauge theory. However, there are 
some major differences, which one encounters in doing 
so. First and foremost, U( 1)2+1 gauge theory is based 
on the principle of local gauge symmetry, absent in 1-D 
quantum mechanics. As a consequence, physical states 
(wave functions of physical particles) have to be gauge 
invariant states. The lattice action S[U] is given by 14511 



S[U] 



flO 



□r 



5elec[f/]+5mag[f/] , 



(23) 



where □( and Dj denote the time-like and space- 
like plaquettes, respectively. The corresponding lattice 
Hamiltonian is given by ll44ll 



^=f^ L lfi + ^Y}l-Re{Uu)]. (24) 



The lattice Hamiltonian has two terms, the electric term 
and a magnetic term. 



lec 



y n. 

2a ^ 'J 

<ij> 

1 



(25) 



The operator /,j counts the number of electric flux 
strings. Its eigenstates are 



lij\hj)^hj\X,j) ,X,7 = 0,±1,±2,. 



(26) 



For each link ij, the states form a complete orthog- 
onal basis. 



I 1^X^1 = 1, 

X=Q.,±l.,... 



(27) 



The magnetic term is built from link operators Uij. It 
has the eigenstates 



Uij\Uii)=Uijpij) 



(28) 



For each link ij, also this basis is a complete orthogonal 
basis. 



J dU\U){U\ = l., 
{U'\U)=b{U' -U) 



(29) 
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The quantum mechanical analogue of the electric flux 
string is momentum P, while the analogue of the link 
is the position X. Sometimes one needs to switch from 
the Bargmann link basis to the electric field string ba- 
sis. In quantum mechanics the scalar product {x\p) can 
be computed from the commutator [X,P] = ; Ti. Like- 
wise, the scalar product (A-|C/) can be obtained from the 
commutator ll20ll 



[l,U] 



-u 



One obtains 



{X\U) = (Uf 



(30) 



(31) 



Secondly, the relation between the transition ampli- 
tude at one hand expressed in terms of the Hamiltonian 
and at the other hand expressed in terms of the path inte- 
gral (in analogy to Eq. (fTTl i). translates in gauge theory 
to 



Mu,.u,„ = mexp[-^T/ri]\Uin 
= / [dU]exp[-S[U]/ri] 



(32) 



Uinfi 



The Hamiltonian requires to choose a gauge and the 
Hamiltonian has been obtained using the temporal 
gauge (Utime-iike = !)■ However, Eq. (|32] | is not true, 
in general. This is due to the fact that a Bargmann 
link state is not a gauge invariant state and hence, 
the amplitude expressed in terms of the Hamiltonian is 
not gauge invariant either However, the amplitude ex- 
pressed in terms of the path integral is gauge invariant, 
because the Haar measure dU is invariant and the action 
S\U] is also gauge invariant. Thus, both amplitudes can 
not be equal, in general. However, one can show that 
they become equal after projection of the Bargmann link 
states onto gauge invariant states ll46ll . that is 



(f/fi|nexp[-i^r]|{/ir 



[dU]txp[-S[U] 



(33) 



Here, the operator fl denotes a projection operator of 
Bargmann states onto gauge invariant Bargmann states, 
which commutes with the Hamiltonian. For a single link 
n is given by 



(34) 



and correspondingly for a multiple link state and 
are group elements). 

Thirdly, we would like to point out that in lattice 
gauge theory, it is customary to identify the measure 
[dU] in Eq. (|32] | with the Haar measure of the group. 
This is not quite right, as there is some normalization 
factor missing. Such normaUzation factor, however. 



cancels out when computing ratios of path integrals as is 
usually done in computing observables via Monte Carlo 
with importance sampling. However, if the path inte- 
gral stands alone, such normalization factor needs to be 
taken into account. In the case of U(l) lattice gauge 
theory, such normaUzation factor has been computed in 
Ref. 1I47I1 . It is given as follows. We consider the Eu- 
clidean transition amplitude for a single time slice aq 
on the lattice under the electric part of the Hamiltonian 
(like in quantum mechanics such normalization factor 
only depends on the kinetic term). The integral of the 
transition amplitude over all final states (which is gauge 
invariant) is normalized to unity. 



dUfi (f/fi I exp [- i?i^iecflo] I f/ir 



[Z«'C/fi]exp[-5elec[f/]|y''''=0° 



+^ dan 
Z — — exp 
2% 



[1 -COs(afi -ttin)] 



= 1 . 

(35) 



Defining A = a/ (;?^ao), and using the Bessel function of 
imaginary argument ll48ll . 

/o(z) = - / t/eexp[zcos(e)] , (36) 

K Jo 

Eq. (|35] | yields 

exp (A) 



HA) 



(37) 



We keep T fixed, and let ^ 00 (which means ao — 
and A which is the continuum limit in time direc- 

tion. However, in space direction we keep a = 1. The 
asymptotic behavior of Z for ao going to zero is given 
by 



^[l-lA-'-^A-^ + OiA-')] 



(38) 



In the limit aq ^ the leading term Z{A) — \J7.%A is 
sufficient to guarantee that the amplitude under the inte- 
gral in Eq. ( |35T l goes over to 5(L'fi — L^in), as should be. 
A numerical simulation ll47ll has shown that in lattice 
gauge theory (with non-zero ao) also the sub-leading 
term is important and can not be neglected. 



IV. STOCHASTIC BASIS 

We want to find a finite set of physically relevant ba- 
sis states {T^ | A* = I7 ■ ■ ■ in Bargmann space. We 
suggest to choose those states randomly, drawn from a 
physically guided distribution. There are several possi- 
bilities in choosing such distribution. For a system with 
a given (Kogut-Susskind) Hamiltonian 9i the physically 
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motivated choice for the distribution is given by the tran- 
sition ampHtude in imaginary time, involving the same 
Hamihonian, 



p(c/) = (f/|exp[:?/r/ri]|f/in) , 



(39) 



where Uia is some suitably chosen fixed spatial lattice 
configuration (Bargmann state). This function is suit- 
able as probability distribution because it is a positive 
function P{U) > 0. The positivity can be seen by ex- 
pressing the transition amplitude in Eq. (|39] l in terms 
of an Euclidean path integral, where the group mea- 
sure is positive and the exponential exp[— 5[C/]/r2] is 
also positive. This choice of distribution has the in- 
convenience of not being analytically computable. Nev- 
ertheless, physically relevant configurations (Bargmann 
states) can be drawn from this distribution by expression 
P{U) as path integral like in quantum mechanics (see 
above) and doing the sampling via Monte Carlo method. 
Still considering the underlying Hamihonian lf{, one 
may consider as alternative the distribution given by the 
transition amplitude from the electric Hamihonian, 



P{U) = {U\exp[-H,iecT/n]\Uir 



(40) 



This distribution analytically is a computable func- 
tion. Moreover, this function is normalized to unity, 
/ dUP{U) = 1. The latter property is due to the fact 
that first the Hamihonian !Ne\ec is a Casimir operator 
for gauge transformation (commutes with all genera- 
tors of local gauge transformations) and second due to 
the left/right invariance of the group measure. The an- 
alytical computation of the distribution P{U) given by 
Eq. ( |40l i is discussed below. Using Monte Carlo tech- 
nique to sample from such distribution P{U), one ob- 
tains the Bargmann states denoted by | T^). Both of the 
above alternatives of distribution involve a time parame- 
ter T, which determines the "width" of the distribution. 
Such time parameter needs to be tuned. As a general 
rule, we used to choose T such that it falls into the scal- 
ing window of eigenvalues (see below). 



V. TRANSITION AMPLITUDES UNDER ELECTRIC 
HAMILTONIAN 

Let us consider the transition amplitude be- 
tween Bargmann link states (for simplicity we con- 
sider first a lattice consisting of a single spatial 
link) under evolution of the electric Hamihonian, 
{U&\exp[—lHeiecT /Ti]\Uin). Its quantum mechanical ana- 
logue is (xfil exp[— i?4in7"/^] l^in)- In elementary quan- 
tum mechanics one learns to compute the latter ampli- 
tude by expanding position states in terms of momen- 
tum states which are eigen states of the kinetic Hamil- 
tonian. This is done via the Fourier expansion theorem. 
One can proceed by analogy to compute the above am- 
phtude in lattice gauge theory. The position x corre- 
sponds to the link variable U. Because U = exp[!a] is a 



periodic function in the interval —% < a < 7t, the conju- 
gate variable of U has discrete values, being the eigen- 
values of electric flux strings, Eq. (l2&t . Thus functions 
of the variable U ~ exp[/a] can be expanded by dis- 
crete Fourier expansion. In case of the Abelian group 
U(l), the Peter- Weyl theorem says that the Bargmann 
link states can be expanded in terms of irreducible repre- 
sentation matrices. In this case, the Peter- Weyl theorem 
is equivalent to Fourier expansion ll49ll . The Peter- Weyl 
theorem holds more generally for groups SU(N) [pOl 
(see Appendix). Applying the Peter- Weyl theorem for 
the group U(l) to the above amplitude yields 

{Un\sxp[~^i,,T/h]\U,n) 

n=0,±l,±2,... ^'^ 

(41) 

In mathematical terms n (running over 0,±1,±2,...) 
denotes the index of the irreducible representation, (t/)" 
also denotes the irreducible representation of group el- 
ement U with representation index (quantum number) 
n. In physical terms, n represents the number of elec- 
tric flux lines (see Eq. (|26]|). The Hamihonian i?4iec is a 
Casimir, which is diagonal in the representation index n. 
Note that the amplitude only depends on the group el- 
ements via the product U^^U[n- By parameterizing the 
link variables via U = exp[!a], one obtains 

<f/fi|exp[-i?4iecr/ri]|f/i„> 



n=0,±l,±2,... 



2a 



n ]cos[n(afi — Oil 



(42) 



A. Construction of gauge invariant states 

The Peter- Weyl theorem can be applied to an arbi- 
trary lattice. It is useful for the construction of gauge 
invariant states. For example, let us take a spatial lattice 
consisting of four links ordered to form a plaquette (see 
Fig.|2] top panel). Then the transition amplitude is given 
by 

(C/fi|exph:?4iecr//^]|f/") 



-54lecr/fi) irrin j fm jfm rrin\ 
l'^127'^237'^43;'^14/ 



n E exp 

17=12,23,43,14 L ;i,-^=0,±l,±2„.. 



cos[n,^(afi-a;;;)] 



-^4 

2a '■' 



(43) 



In order to make the amplitude gauge invariant, we carry 
out the group integral / d'^^i — ^ J^^d^i at the nodes 
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1,2,3,4, as 

(c/fi|n /n]\u 
1 

2k 



1 • • • 



271 



X n L ^'^p 

i7=12,23,43,14 «,-^=0,±l,±2,... 

X cos[«,Xafi-(a;^ + p,-P^))] 



^'^^7- 2 
2fl 



(44) 



The group integral introduces Kronecker delta functions 
at each vertex. E.g., at vertex j, one obtains dnjj.nji,- The 
number of ingoing flux lines equals to the number of 
outgoing flux lines. As we consider here the absence of 
any charge, this rule represents Gauss' law. By defining 
the plaquette angle 



6plaq = 0Cl2 + + a43 + "14 
A9piaq = 9piaq - 9p"aq , 



(45) 



we obtain the final expression of the gauge invariant am- 
plitude, 

(u^jcxp[-H.i,,T/ri]\uZ) 



E 

«=0,±1,±2,... 











exp 


2a 


cos 


nAQpiaq 



(46) 



Here n denotes the number of closed plaquette loops on 
top of each other The result is built from plaquettes 
which are closed loops of consecutive link variables, 
forming the smallest non-local gauge invariant objects 
on the lattice. Note: Opposite signs of n correspond to 
plaquettes of opposite orientation. The eigenvalue of the 
electric field corresponds to the contribution from « 
plaquette loops. It is also important to note that the re- 
sult only depends on the number of plaquette loops and 
the difference between initial and final plaquette angles. 

This result can be generalized to any 2-D spatial lat- 
tice including A^p plaquettes, as following: 



{Ut^\exp[^iH^i^,T/ri]\ul 

Np f 

n E exp 

P=l I- np=0,±l,±2,... 



X COS rip.AQp 



2a sraph 



(47) 



For example, in the case of a spatial lattice composed 
of four plaquettes (Fig. [2] middle and bottom panels) the 
gauge invariant amplitude becomes 



{Ul'jexp[-!}ii,,T/n]\ut 



X cos 



E 

ni ,n2-^3-«4=0-il-i2 
plaq 



exp 



"1^91254 - 



«2 AeP3g''5 



2a ^S'^'ph 

-«3 AePggV 



- «4 A9; 



plaq 
'4587 



(48) 



U43 

4, ,i£ ,3 



U,4 



U23 







7, 








■ U45 


■^8 ■ 
5 




, ^5 


1^6 


^4. 



1^9 



1^6 



^2 

9 







"3 






' '"3 nj^ 


n4 




^3 


J4 




nj 5 


"2 




>'ni 


'"2 








— 1 1 



1^3 
7, 



6 <^ 4* 



8 



"3 



"1 



04 



"2 



1 



FIG. 2: Spatial lattice. Top panel: The state of a spatial pla- 
quette described by Bargmann link states \ Uij). Middle panel: 
the same scheme for a lattice including four plaquettes. Bot- 
tom panel: The state of lattice can also be expressed in terms 
of electric flux states, by assigning a link flux number nij and a 
direction of flux lines between neighbor sites i and j. A gauge 
invariant state obeys Gauss' law at each vertex, i.e., the total 
number of incoming minus outgoing flux lines is zero. On a 
2-D lattice, a gauge invariant state can be characterized by as- 
signing a plaquette flux number iipi plus an orientation to each 
plaquette. 



where 



4raph = 2«? + 2«i + 2«^ + 2«2 + 

+ («! - nif' + (ni - n^f + (n2 - n^f + (na - n^f- . 

(49) 



^graph denotes the eigenvalue of the electric field op- 
erator of the graph constructed from plaquettes fill- 
ing the lattice, each with counter-clockwise orienta- 
tion. Analytic results have been obtained also in three- 
dimensional spatial lattices, where instead of plaquettes 
one has cubic states. Because we present numerical re- 
sults in 2-D lattices, we do not give details of 3-D here. 
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TABLE I: Electric Hamiltonian. Comparison of eigenvalues 
of effective Hamiltonian versus exact eigenvalues. Lattice 2^, 
spacing a = aQ = I, time parameter in distribution P{U) is T = 



0. 1, coupling g = 1, dimension of stochastic basis A'ba; 



■■32. 









p (exact) 


Rel. error 





1.0000000000 


0.0000000 


0.0 








1 


0.8187308431 


1.9999989 


2.0 


5.5x 


10^ 


-7 


2 


0.8187306875 


2.0000008 


2.0 


4.0x 


10^ 


-7 


3 


0.4493290180 


7.9999988 


8.0 


1.5x 


10- 


-7 


4 


0.4493289101 


8.0000012 


8.0 


1.5x 


10^ 


-7 


5 


0.1652988898 


17.9999999 


18.0 


5.5x 


10- 


-9 


6 


0.1652988667 


1 a f\f\f\r\f\ 1 o 

18.0000013 


18.0 


7.2x 


10- 


-8 


7 


0.0407622060 


31.9999995 


32.0 


1.5x 


10- 




Q 

o 


U.U4U /DZ1V41 


jZ.UUUuUZ4 




7.5 X 


10- 


-8 


9 


0.0067379553 


49.9999876 


50.0 


2.5 X 


10- 


-7 


10 


0.0067379496 


49.9999961 


50.0 


7.8x 


10- 


-8 


11 


0.0007466034 


71.9997633 


72.0 


3.3x 


10- 


-6 


12 


0.0007466023 


71.9997789 


72.0 


3.1x 


10- 


-6 


13 


0.0000554676 


97.9971144 


98.0 


2.9x 


10- 


-5 


14 


0.0000554663 


97.9973409 


98.0 


2.7 X 


10- 


-5 


15 


0.0000027819 


127.9235530 


128.0 


6.0x 


10- 


-4 


16 


0.0000027792 


127.9332884 


128.0 


5.2x 


10- 


-4 


17 


0.0000001401 


157.8084005 


162.0 


2.6 X 


10- 


-2 


18 


0.0000001134 


159.9218144 


162.0 


1.3x 


10- 


-2 


19 


0.0000000873 


162.5307945 


200.0 


1.9x 


10- 


-1 


20 


0.0000000810 


163.2828884 


200.0 


1.8x 


10- 


-1 



VI. NUMERICAL RESULTS 




A. Test of effective electric Hamiltonian 

In order to test the effective Hamiltonian, first we 
consider the electric Hamiltonian. This is a good test 
bed, because the physics of the electric Hamiltonian can 
be computed analytically. First and as an example, for 
the case of a 2^ lattice, we present the eigenvalues !D„ 
of the matrix M and the energy eigenvalues E„ of the 
effective Hamiltonian versus the exact eigenvalues, as 
well as the relative errors (see Tab. The following 
behavior is observed. The ground state, corresponding 
to «piaq = (zero plaquette loops), has no degeneracy. 
All excited states are doubly degenerated, correspond- 
ing to «piaq = ±1,±2, The relative error stays in 

the order of 10^^ for the first 5 levels, and then grows 
exponentially with increasing level number n. The rel- 
ative error reaches the order of 1 for « = 20. Beyond 
that the spectrum of the effective Hamiltonian drowns 
in numerical noise. This means that there is an energy 
window, going from £ = to £ = 165 where the effec- 
tive Hamiltonian gives meaningful results. Such energy 
window is correlated with the behavior of diagonal ele- 
ments T>„ of the transition matrix M. The computations 
have been done with double precision, i.e., 15-16 digits. 
One observes that the upper bound of the energy win- 
dow is reached, when the value of 'Diq w 0.81 x 10^^ 
has 8 significant digits being exactly half of the 16 dig- 



FIG. 3: Electric Hamiltonian. Eigenvalues of effective Hamil- 
tonian are compared to exact eigenvalues (0, 4.5, 6.75, 9.0, 
1 1 .25, 13.5 and 18.0). Top panel: spatial lattice 3"^, imaginary 
time P = 2.0, spacing a= I, coupling constant g = 1.5, asym- 
metric factor ^ = 20, number of stochastic basis A'ba.sis = 400. 
Middle panel: lattice 8 x 8, P = 4.2, a = I, g = 1.5, ^ = 5, 
Miasis = 1000. Bottom panel: The same results as the mid- 
dle panel along with relative errors. Notice that the extension 
of horizontal parts in the top and middle panels indicates the 
eigenvalue degeneracy, which grows with the lattice size. 

its of internal arithmetic precision of the computer The 
existence of such energy window has been observed in 
all cases. The size of the energy window depends on the 
following parameters: (i) the internal numerical preci- 
sion of the computer, and (ii) physical parameters, like 
coupling constant g, lattice size a, and transition time T. 
Notice that the transition time is equivalent to p, which 
is the inverse of temperature (in all figures we use P in- 
stead of T). The internal numerical precision can be 
viewed as a scale of experimental observation, similar to 
the wave length of light in a microscope. The obtained 
results depend on the internal scale of the experimental 
apparatus, which in our case is the internal arithmetic 
precision of the computer If we want to increase the 
size of the energy window we can (i) use higher numer- 
ical precision, which may be computationally costly, or 
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FIG. 4: Electric Hamiltonian. Scaling windows of energy 
eigenvalues of effective Hamiltonian are presented (the full 
line is a guide to the eye). The regime where the line is hor- 
izontal (energy scales in P with exponent zero) represents a 
scaling window. Top panel: Spatial lattice 2^, a = 1, ^ = 20, 
g = 1.0, A'basis = 200. Middle panel: Spatial lattice 3^, a = 1.0, 
^ = 20, g = 1.0, A'basis = 200. Bottom Panel: Spatial lattice 
8^, a = 1, ^ = 5, g = 1.5, A'basis = 1000. 



(ii) choose physical parameters appropriately. 

One finds that larger lattices lead to an increased de- 
gree of degeneracy in the energy spectrum. Such de- 
generacy comes from small closed loops, which located 
anywhere on the lattice, will give the same electiic en- 
ergy. Such high degeneracy will be lifted when taking 
the magnetic term into account. With increase of spatial 
lattice size an increase of the stochastic basis dimen- 
sion is required. Also the tuning of time parameter T in 
distribution P{U), (Eq. |40l i makes that T changes with 
lattice size. 

The above results all correspond to a particular value 
of transition time T in the transition amplitude. How- 
ever, the exact physical energy levels must be indepen- 




0.05 



FIG. 5: Scaling of wave functions for the electric Hamilto- 
nian. Expansion coefficients of wave functions comsponding 
to eigenvalues in Fig. |4] (top panel) in terms of first stochastic 
basis. Spatial lattice 2 x 2, a = 1, ^ = 20, g = 1, A'basis = 200. 



dent of the transition time T . The dependence of energy 
levels under variation of the parameter P(= T) has been 
investigated. Results are shown in Fig. |4] for three spa- 
tial lattices. The top panel (the 2^ lattice) shows that 
the ground state is a flat line as a function of p, in a 
wide range of p (window from p = 0.2 up to p = 5), 
meaning that the ground state energy is independent of 
the transition time p. In other words, the ground state 
energy scales like p^ with exponent y = in a finite 
time window. A similar scaling window is observed also 
for the higher excited states. However, with higher en- 
ergy, the size of such scaling window becomes smaller. 
For energy E > 100, the scaling window is no longer 
visible. We found that the size of the scaling window 
S„ can be described approximately by an exponential 
law S„ °^ exp[—aE„]. Such behavior of decreasing scal- 
ing windows can be understood from the property that 
e\p[—^HT/H] projects onto the ground state for large T 
(Feynman-Kac theorem). Higher levels become expo- 
nentially suppressed by the dominant ground state and 
can survive only for shoit times T. The middle and the 
bottom panels of Fig.|4]also display the scalling window 
of 3^ and 8^ spatial lattices. One observes that scal- 
ing windows are observed in a spectrum with higher de- 
generacy. For the 8^ lattice the distribution P{U) is ex- 
pressed via the path integral, while for the 2^ and 3^ lat- 
tices it is given by the analytic expression (see Eq.l47]i. 

The scaling behavior observed in the energy eigen- 
values should be manifest also the corresponding wave 
functions. In particular, we have studied (e^|<I>„), i.e. 
the expansion coefficient of wave function <t>„ in terms 
of the stochastic basis function e^,. As two examples, 
such expansion coefficients, expanded in terms of the 
first basis function, are shown for the lattice sizes 2^ 
and 8^ in Fig. |5] and Fig. |6l respectively. Indeed, one 
finds scaling behavior also in such expansion coeffi- 
cients. Like the size of the energy scaling window de- 
creases with increasing energy E„, also the wave func- 
tion has a scaling window, and its size also decreases 
with the level index n of energy. Moreover, the eigen 
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FIG. 6: Electric Hamiltonian. Scaling window of the expan- 
sion coefficients of wave functions, corresponding to eigenval- 
ues in the bottom panel of Fig.|4] Spatial lattice 8 x 8, a = 1, 
^ = 5, g = 1.5, A'basis = 1000. Notice that yti = 1 in the top 
panel relates to the vacuum state, and k2, k^, in the other 
panels are indices of the second, the third and the forth ex- 
cited levels. Due to the high rate of degeneracy of these states, 
only some of them have been shown here. 



function degeneracy is observed for the 8^ lattice, and it 
is increased with the lattice size. 

Thermodynamics. The quality of the energy spec- 
trum of the effective Hamiltonian can also be seen from 
the thermodynamical functions. In the case of the elec- 
tric Hamiltonian, the energy spectrum Eo,Ei can be 
computed analytically. Thermodynamical functions can 
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FIG. 7: Thermodynamical functions computed from exact 
eigenvalues (full line) versus those from eigenvalues of ef- 
fective electric Hamiltonian for a 6 x 6 lattice with a = I, 
^ = 5 (ao = 0.2), g = 1.0 and Wbasis = 600. Top panel: Aver- 
age energy t/(P). Bottom panel: Specific heat C(p). In lower 
temperatures (higher values of p) there is a reasonable agree- 
ment. 



be expressed via those energies. The partition function 
Z is denoted by 

Z(p) = Tr [exp(-pJ^)] = ^exp [-p£„] , (50) 
and the average energy t/(|3), e.g., is given by 

aiogz(p) 1 



um = - 



^£„exp[-p£„] . (51) 



Likewise, one can compute free energy, entropy and 
specific heat. For example. Fig. [T] shows the results of 
the average energy C/(p) and the specific heat C(p) for 
the electric Hamiltonian of a 6^ lattice. In this figure 
we have compared the results from the effective Hamil- 
tonian with the exact Hamiltonian. In general one ob- 
serves that agreement is good in the regime of large P, 
i.e., the low temperature regime. For small values of p 
some disagreement becomes visible, reflecting the fact 
that the precision of higher energy levels of the effec- 
tive Hamiltonian is limited (their scaling windows go to 
zero). 



VII. RANDOM PERTURBATION OF ELECTRIC 
HAMILTONIAN 

Below we will present results for the effective Hamil- 
tonian of the full Hamiltonian. In order to better un- 
derstand the propagation of numerical errors, we have 
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considered the electric Hamiltonian plus a small random 
perturbation, which shall simulate numerical errors. We 
have then computed the spectrum of the corresponding 
effective Hamilton, and compared the cases with and 
without the random perturbation. The results are shown 
in Fig. [8] One clearly observes that the random pertur- 
bation reduces the number of states which show scaling 
and also reduces the size of the scaling windows. 
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FIG. 8: Electric Hamiltonian. Influence of artificial random 
relative error 5 in transition matrix on scaling behavior of en- 
ergies. A 2^ spatial lattice. Top panel: a = 1, ^ = 20 (lattice 
asymmetry parameter), g = 1, A'basi.s = 200, 5 = 0.01 (no er- 
ror). Bottom panel: the same graph with error (S = 0.01). 



VIII. TRANSITION AMPLITUDES UNDER FULL 
HAMILTONIAN 



ant projected states is written as 

M,,v(r) = (t/^|n exp[-:?4iirA]|f/v) 



J[dU] exp[+S, 



mag 



[u]/n]eM-s[u]/h]\''^l'l 



(52) 



Here Uf, denotes the stochastic basis of Bargmann states 
on the lattice (in analogy to the quantum mechanical 
stochastic position states X/j). We have drawn these 
states from the distribution corresponding to the elec- 
tric Hamiltonian, given by Eq. ( l40t . We have preferred 
this choice because this function is analytically com- 
putable, which means a numerical effort being substan- 
tially smaller than that for the distribution involving the 
full Hamiltonian, Eq. (3% . We construct a matrix of 
transition elements between normalized stochastic basis 
states. We proceed in analogy to quantum mechanics 
(Eq. [TtI i. Finally, we diagonalize such matrix and ex- 
tract eigenvalues and wave functions (Eqs.[T8]|2T]i- 



25 
20 
15 
10 
5 








8x 8 Lattice 








a = 1.0 








5 = 20 








g = 2.5 

. =600 

basis 






1 


2 


P 3 4 


5 



o 



0.2 




8x 8 Lattice 


0.15 




a = 1.0 
5 = 20 


0.1 




g = 2.5 
k . =600 
& basis A 


0.05 


















Let us consider the gauge invariant transition ampli- 
tude under the full Hamiltonian (in analogy to quan- 
tum mechanics (see Eq.lfTTI)). Although this can be ex- 
pressed in terms of a path integral with the lattice ac- 
tion (Wilson action), this is numerically not suitable, 
because Monte Carlo with importance sampling only al- 
lows to compute ratios of transition amplitudes. Hence, 
like in quantum mechanics (Eq.(fT4li) we factorize the 
above amplitude into two terms, one being analytically 
computable and the other one being given by the ratio 
of transition amplitudes computable via the method of 
Monte Carlo. Therefore, the transition amplitude under 
evolution of the full Hamiltonian, between gauge invari- 



FIG. 9: Full Hamiltonian: Scaling window of the spatial 8x8 
lattice with a = I, g = 2.5, ^ = 20, A'basis = 600. To panel: 
ground state energy. Bottom panel: expansion coefficient of 
the ground state wave function. 

Full Hamiltonian: influence of the magnetic term. 
Taking into account the magnetic term generates the ef- 
fective full Hamiltonian. First results on scaling of en- 
ergy eigenvalues and wave functions of low-lying states 
are shown in Fig. |9] These results correspond to a 
8^ lattice. The results show scaling windows for the 
ground state energy (top) and the expansion coefficient 
of the ground state wave function (bottom). Comparing 
with the scaling observed in the electric Hamiltonian, 
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only one level shows scaling, and the scaling window 
is smaller. This can be understood from the fact that the 
ratio of matrix elements in Eq. (l52t has been determined 
via Monte Carlo from path integrals, which carries sta- 
tistical errors in the order of a few percent. We have 
seen that errors of such order of magnitude de reduce the 
scaling behavior (see Fig.[8]l. From this observation we 
conclude that the numerical resolution of energy levels 
and scaling windows in the full Hamiltonian is essen- 
tially determined by the statistical error occurring in the 
numerical calculation of the ratio of matrix elements. 
Results with better statistics and a larger stochastic ba- 
sis are being called for 

IX. DISCUSSION 

The construction of a Hamiltonian in lattice gauge 
theory faces the following problems; 

(a) Although there exists the lattice Hamiltonian, this 
alone has not proven to be useful in computing viable 
results for physical observables. We suggest here that 
the best one can do is to construct an effective Hamil- 
tonian. This is meant to be a Hamiltonian which de- 
scribes physics in a finite window (e.g. a window of 
low energy). This has some analogy to the idea of of the 
Wilson-Kadanoff renormalisation group, where a renor- 
malized Hamiltonian is constructed, which is valid at 
some critical point, but distant from the critical point 
has no physical meaning. 

(b) Conventional Hamiltonian methods used in 
atomic physics, nuclear physics, condensed matter and 
particle physics consider the given Hamiltonian J-C and 
compute matrix elements ((]),■ from such Hamil- 
tonian. In contrast, here we consider a function of the 
Hamiltonian exp[— i^/^r//;]. This has the following ad- 
vantages: First, ^ is mathematically a more singu- 
lar object than exp[— i?/]. This can be seen at hand 
of the simple quantum mechanical example of the ki- 
netic Hamiltonian i?4in- The matrix element (yli^nl^i^) 
is a derivative of a 5-function, while the matrix element 
(yl exp[— i?4in] is a smooth, differentiable and rapidly 
falling-off function. This means, the non-linear expo- 
nential function smoothes out singularities of an oper- 
ator Second, matrix elements of the exponential func- 
tion exp[— i^/r//;] can be evaluated by using the path 
integral. Third, and mostly important, in contrast to J-[, 
the operator function exp[— i?/^r/rj] contains a parame- 
ter r, which is redundant for the physical spectrum, i.e. 
using any value of T one should obtain the same spec- 



trum of In numerical simulations such redundancy 
gives additional information about errors in the follow- 
ing way: The results do depend on the value of T . They 
depend also on physical parameters, such as coupling, 
lattice size and lattice spacing etc. They further depend 
on approximation parameters, like the number of equi- 
librium configurations used in the path integral, the size 
of the stochastic basis etc. Last but not least they depend 
on the internal precision used in the computer Here 
we can turn the dependence on T into an advantage: (i) 
The T-dependence of the energy spectrum (or better of a 
number of low lying energy eigenvalues) can serve as a 
measure of error of the calculation. In the best case the 
energy eigenvalues become T independent. This hap- 
pens in the so-called scaling windows. In the worst case 
they are strongly T-dependent, meaning that these re- 
sults are unphysical. (ii) We can use the scaling window 
to tune the time parameter occurring in the distribution 
P{U), which generates the stochastic basis. 

(c) Stochastic basis. In our opinion, the construction 
of suitable basis is the most important step in order to 
compute physics from a Hamiltonian. Such basis is built 
on two principles: first a random pick and second, a 
physical principle to guide the search. The so-called 
stochastic basis is built in close analogy to the equi- 
librium path configurations computed via Monte Carlo 
importance sampling to solve Lagrangian lattice gauge 
theory. 

(d) In contrast to Lagrangian lattice gauge theory, 
where gauge symmetry is manifestly conserved in the 
path integral via the group measure and the (Wilson) 
action, in the Hamilton formulation gauge invariance 
of states and amplitudes has to be imposed (via Gauss' 
law). Technically, much work is required to construct 
such gauge invariant states. One expands the link 
Bargmann states into irreducible representations using 
the Peter- Weyl theorem. Gauge invariant states are then 
constructed by doing the group integral of local gauge 
transformations (at each vertex). As a result one en- 
forces Gauss' law at each vertex. Here, we have shown 
how this can be done in the case of U(l) gauge theory. 
This can be generalized to non-Abelian gauge symme- 
try. For example, Burgio et al. lIsUl have shown how to 
construct a gauge invariant Hilbert space for the gauge 
group SU(2). 
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X. APPENDIX A : PETER- WE YL THEOREM FOR 

SU(N) 

We consider the Lie group SU(N). Irreducible repre- 
sentations play a crucial role in the Peter- Weyl theorem 
in the case of U(l) and also for SU(2). This is also true 
in general for SU(N). Let us consider irreducible repre- 
sentations in the notation of Young diagrams character- 
ized by a partition and denoted by 

{v} = {vi;...;VAf_i} , where vi > ••• > Viv_i . (53) 

It is convenient to define the number 

V = 1^ V, . (54) 

1=1 

Elements of the group manifold can be conveniently 
parametrized by 



: exp 



(55) 



where Ti , . . . denote the group generators in the 

fundamental representation. A group element in the 



irreducible representation characterized by index (set) 
{v} is given by a matrix 

Dllhm, (56) 

where a,b aie the matrix indices running from 1 to the 
dimension |v| of the irreducible representation. 

In order to compute the transition amplitude of the 
electric Hamiltonian, one can use the following corol- 
lary of the Peter-Weyl theorem. Let C denote the 
quadratic Casimir operator, 

C= ^ t]. (57) 

Let Cv denote the eigenvalue of C in irreducible repre- 
sentation {v}. Then as a corollary of the Peter-Weyl 
theorem, one can express a matrix element of an opera- 
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tor function /(C) by 

{v}a.b=i 

^ |v|/(Cv)rr[DM(j^-i^^')]^ 

V 

L |v|/(Cv)x^'^(^^-'?^')- (58) 

V 

The sum over {v} runs over all partitions {v} = 
{vi;...;VAr_i} withvi > ••• > Vat-i and 

xW(^) = rr[DM(^)] (59) 
denotes the group character 



4x 4 Lattice 
a= 1.0 




4x 4 Lattice 
a = 1.0 
^ = 20 
g = 2.5 



N. . =200 




